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Abstract 



We discuss the universal dynamics of elastic interfaces in quenched random me- 
C^ , dia. We focus in the relation between the rough geometry and collective transport 

properties in driven steady-states. Specially devised numerical algorithms allow us 

to analyze the equilibrium, creep, and depinning regimes of motion in minimal 
^ ', models. The relevance of our results for understanding domain wall experiments is 
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Resume 

Approches Numeriques pour les interfaces elastiques en milieu aleatoire 

Nous discutons les principaux resultats obtenus sur les proprietes universelles de la 
dynamique des interfaces elastiques en milieu aleatoire. Une attention particuliere 
^T ■ sera dediee a la relation entre la geometric rugueuse de I'interface en mouvement 

^r^ . et ses propietes de transport collectif. Les approches numeriques developpees per- 

mettent de decrire les proprietes d'equilibre, la dynamique de creep et la transition 
de depiegeage de I'interface. Nous discutons aussi la pertinence de nos resultats dans 
ISl • les experiences sur la dynamique des parois. 
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1 Introduction 



The out of equilibrium interplay between quenched disorder and elasticity in 
uniformly driven interfaces is at the root of the universal dynamical response 
displayed by very diverse physical systems. Examples are magnetic [1,2,3,4,5] 
or ferroelectric domain walls [6,7,8,9,10], contact lines in wetting [11,12], frac- 
tures [13,14], and earthquakes [15]. In this paper we discuss the basic phe- 
nomenology that emerges by solving, with specially devised numerical meth- 
ods, the minimal models proposed for capturing such dynamical behavior. 

Universal dynamical properties can be captured, both qualitatively and quan- 
titatively, by rather simple models. To be concrete we will focus on the paradig- 
matic driven quenched-Edwards- Wilkinson (QEW) universality class, mini- 
mally described by 



7(9(m(x, t) = cdlu{x,t) + Fp{u,x) + f + ^(a;, t). (1) 

This equation models the overdamped dynamics of an elastic interface with 
a univalued scalar displacement field u{x,t), with x a vector of dimension d, 
such that the interface is embedded in a space of dimension D = d + 1. The 
elastic approximation made above, assumes small deformations such that the 
elastic part of the energy is well described by Hei = (c/2) /^ dx'^[Vu]'^. We are 
thus ignoring "plastic" deformations such as overhangs or pinched-ofF loops 
that might appear in real interfaces, and assuming that elastic interactions are 
short-ranged, with a stiffness constant c. The contact with a thermal bath at 
temperature T is modeled by a Langevin noise, such that {ri{x,t)ri{x',t')) = 
2kBT'~f6{x — x')S{t — t'). Finally, the interface is coupled to a uniform driving 
force / and to a quenched pinning force Fp, which arises from the disorder in 
the host materials. We will consider a non-biased pinning force characterized 
by its disorder-averaged correlator 



Fp{u, x)Fp{u', x') = A(m - u')6{x - x'), (2) 

with A(x) a short-ranged function, of range Vf (determined by the domain 
wall width in real interfaces with point disorder). For the so-called random 
bond (RB) case, the elastic line moves in a random short-range correlated 
potential. The corresponding pinning force is Fp{u,x) = —duU{u,x), where 
U{u, x) is the random potential, and thus /„ A{u, x) = 0. For the random field 
(RF) case, U{u, x) is a random walk as a function of u, with diffusion constant 
/„A(m,x)>0[16]. 

The model just defined is minimal, and in order to compare with experi- 
ments other ingredients might be considered. For example in charge den- 
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Figure 1. Linking transport and geometry, (a) Snapshot of a domain wall in a 
two dimensional ferromagnet [27]. (b) Typical velocity-force characteristics, (c) 
Crossover lengths £opt and £av representing the optimal excitation and the deter- 
ministic avalanches, respectively, (d) Geometric crossover diagram. 

sity waves and vortices, the elastic structure is periodic [17,18], while in 
fracture [14,19,20,21,22] and wetting [23,24] the elastic interactions are long 
ranged. Moreover, anharmonic corrections to the elasticity or anisotropics can 
be also relevant [25,26]. Remarkably, all these different universality classes 
(with different critical exponents) share the same basic physics of the QEW 
model discussed in the following sections. 



2 Transport and Geometry of Driven Interfaces in Random Media 



Disorder makes the interface dynamics very rich. On one side, the interface 
appears rough, both in absence and in presence of drive. On the other side, 
the response of the system to a finite drive is strongly nonlinear and, at zero 
temperature, motion exists only above a finite threshold, the so-called critical 
force fc- As summarized in Fig. 1 collective dynamics is observed in different 
regimes: just above fc the motion is very jerky, and displays collective rear- 
rangements or avalanches of a typical size £av Below fc motion is possible only 
by thermal activation over energy barriers and the interface slowly fluctuates 
backward and forward. This "futile" and reversible dynamics takes place up 



to a characteristic sizeLJ iopt above which energy barriers do not exist before 
the next lower-energy metastable state and the interface moves only forward 
pulled by the finite drive /. As / — )■ both i^pt and the energy barriers diverge. 

2.1 Three reference stationary states 



Universality manifests, both in the transport and in the geometry of the mov- 
ing interface, through the existence of robust critical exponents. These describe 
the rate of power-law divergences of important quantities as the control pa- 
rameter / approaches three special states: (i) the equilibrium (/ = 0); (ii) 
the depinning (/ = /c, T = 0), and the fast-flow (/ ^ fc). Beyond some 
microscopic length-scale [28,29] these three steady-states have the pecuharity 
that the geometry is self-affine with their characteristic roughness exponent: 
(i) Ceq, (ii) Cdep, and (iii) Cff- Self-affinity means that length x and displace- 
ment u can be rescaled as x' = ax and u' = a~^u into a new interface u'{x') 
which is statistically equivalent to u{x), namely u{ax) ~ a'^u{x). Transport 
properties in these three points are very different: At equilibrium, the mean 
velocity is zero and the dynamics is glassy: at small length scales, we observe 
a fast futile motion, while in order to observe a rearrangement on a large scale 
size / we need to overcome a barrier -£";,(/) growing as E^ r^ l^, with 6 a positive 
and universal exponent. Dimensional analysis suggests that 6 = d — 2 + 2(^eq, 
which is confirmed by accurate numerical studies [30]. At the zero tempera- 
ture depinning transition the velocity vanishes as v{f,T = 0) ~ (/ — /c)^ for 
f > fc while V = for / < fc- At finite temperature, this sharp transition 
is rounded and the velocity behaves as v{fc,T) ~ T^. This is the so-called 
thermal rounding regime. At large force, / S> fc, in the fast-flow regime, we 
recover the linear response f ~ /. Here impurities generate an effective ther- 
mal noise on the interface with Tgff — T ~ A(0)/f (the disorder strength A(0) 
is defined in Eq.(2)). Therefore, the fast-fiow roughness corresponds to the 
Edwards- Wilkinson roughnessLlI Cff = (2 — d)/2. The three reference state are 
schematically represented in Fig. 1(b). 



2.2 Connecting the three reference stationary states: creep and depinning 



How does the interface behaves in between (i.e., < f < fc and / > fc) these 
three reference steady-states? Let us consider some instantaneous snapshots 



1. The subscript opt stays from "optimal" because ^opt is the size of the jump 
associated with the optimal barrier that the interface should overcome in order to 
find a new metastable state with a smaller energy. 

2. In real experiments we expect non-linear terms to become relevant and change 
the Edwards- Wilkinson into the Kardar-Parisi-Zhang universality class 



of an interface as obtained from an experiment or a numerical simulation. 
Fig. 1(a) is an illustrative experimental example for a ferromagnet [27]. Let us 
assume that the longitudinal size of the snapshot is large enough to contain 
all the relevant characteristic length-scales, and that the driven interface is 
already in a stationary regime so that the memory of the initial condition is 
lost. By analyzing the snapshot, what can we say about the interface state of 
motion? 

One could imagine a naive scenario where the dynamic roughness exponent 
( varies continously upon increasing the driving force, from its value at the 
equilibrium ^eq, to ^dep, and finally to Cff- Actually, this is not true, and the 
interface geometry can be instead described by only these three exponents and 
by the two crossover lengths £opt and £av The crossover or dynamical phase 
diagram is schematically shown in Fig. 1(d). At small temperatures and be- 
low threshold, < f < fC the interface is in the ultra-slow creep regime. At 
small length-scales i < iopt the interface looks like an equilibrated interface 
(with an exponent Ceq)- For intermediate scales, iopt < i < ^av, we expect 
the same roughness exponents as those of the depinning transition, (^^ep- Fi- 
nally, at the largest length-scales, £av < ^, we expect to measure the fast-flow 
roughness exponent C,±i. Let us observe that the large length scales are con- 
trolled by non-equilibrium roughness exponents. This shows that even for very 
small forces /, the system is far from equilibrium [16,31,32], in contrast with 
the initial physical pictures that assumed metastable configurations indistin- 
guishable from the equilibrium ones. Analogously, for / > /c we observe that 
at short length-scales i < £av the roughness exponent is Cdep, while for i > £av 
one has Cff- 

The dynamical phase diagram of Fig. 1(d) thus allows us to get important 
information from snapshots such as the one of Fig. 1(a): (i) It can tell us first if 
the interface of the snapshot was in the f < fc creep, in the / ^ /c depinning, 
or in the f ^ fc fast-flow regimes; (ii) The actual values of the different 
roughness exponents (obtained by fltting, for instance, the structure factor 
Sq) guide us in the search for a universality class; (iii) Looking at snapshots 
for different forces could give access to the critical behavior of £opt and fav, 
namely i^pt ~ f^"""^ and ia_-v{T = 0) ~ (/ - /c)"''''"^ At flnite temperature, £av 
remains flnite below threshold but quickly diverges as / — ?> [16]. 

Interestingly, both £av adn £opt have a "double" physical meaning: besides be- 
ing roughness crossover lengths, they control the non-linear collective trans- 
port properties of the interface. Indeed, in the depinning regime / > fc, 
T = 0, the jerky motion can be characterized by a velocity autocorrela- 
tion with flnite spatial and time ranges £av and i^v^"^"'^, respectively (with 
Ziep the dynamical exponent). Since the width of the avalanches is expected 
to grow as ^av^"""" the velocity can be thus estimated as v ~ ^av^"^''"^"^''- On 
the other hand f ~ (/ — fc)^ ~ iav~ '^''^, yielding the hyperscaling relation 



P — '-'depl^dep ~ Cdepj- 

In the f < fc creep regime, three different dynamical steps can be iso- 
lated (see Sec. 4. 2): (i) Starting from a deep metastable state, the interface 
explores the neighborhood jumping back and forth (futile dynamics); (ii) 
This continues until a saddle configuration is found. The saddle configura- 
tion will differ from the initial metastable state on a typical length i^pt', (iii) 
Finally, from the saddle configuration the interface relax deterministically to 
a new and deeper metastable configuration. The size of the optimal excita- 
tion iopt can be estimated by balancing the gain in energy of being pinned 
in a deep metastable state, -Emetast(^)! with the gain in energy of moving the 
interface forward ~ / ('Usaddie(^) "~ ""metast (^) ) • At equilibrium we have that 

Wsaddle(^) - Mmetast(^) ~ ('^"'^ and that £;pinned(^) ~ ^^ with 9 = d-2 + 2Ceq. 

The balance gives ^opt ~ f~'^'"^, with z/gg = 1/(2 — ^eq), which means that £opt 
diverges when / — )■ 0, as expected. Moreover, numerical simulations give a 
clear evidence that the energy landscape is characterized by a unique energy 
scale, and that the energy difference between neighbor metastable states is 
equal to the energy barrier separating them [30]. Thus, we expect the barrier 
that the interface has to overcome to escape from a deep metastable state to 
grow as £opt • Using the Arrhenius activation, we recover the creep formula: 

v^exp{-Cf-^^^/T) (3) 



where /i = (2Ceq — l)/(2 — ^gq) is an equilibrium exponent. This formula has 
largely been verified by experiments on magnetic domain walls [1,4,5]. When 
the barriers are too small compared with T, the velocity is no longer described 
by the Arrhenius law, leading to a thermal rounding behaviour v ~ £s.-v~ '^^^, 
with £av ~ j'-^'^dep/i3 g^^ j- — j^^ Yox harmonic short range elasticity it can be 
shown that z/dep = 1/(2 — Cdep) (Statistical Tilt Symmetry relation [33]). 

As summarized in Fig. 1, geometry and transport are thus closely related. 
In the following sections we explain how the above phenomenology can be 
numerically obtained. 



3 Method and Observables 



The main difficulty with Eq. (1) is the non-linearity introduced by the dis- 
order, which breaks the translational symmetry. Indeed, in the absence of 
disorder Eq. (1) becomes the Edwards- Wilkinson equation which is exactly 
solvable [34,35]. With disorder, the mean field model, valid above the upper 
critical dimension rf^c = 4 was solved by Fisher [36]. Advanced analytical tech- 
niques such as the functional renormalization group (FRG) allows to obtain 



expansions below but close to d^c [16]. Numerical approaches thus appear as 
a valuable and necessary theoretical tool. From now on we will focus on the 
d = 1 QEW model, which is: (i) experimentally relevant for interfaces in thin 
films, (ii) a stringiest case for testing analytical approaches in principle valid 
only close to d-ac (iii) simpler to tackle numerically at large length scales. 

Numerically, it is convenient to discretize the interface in the x-direction, 
keeping u{x,t) as a continuum variable! ^ I. The center of mass velocity for an 
interface of size L is defined as, 



1 ^-^ 
vit) = -J29Mx,t), (4) 

^ x=0 



that, given Eq. 1 and r^ = c = 1, is nothing else but the spatial average of 
the instantaneous total forces acting on the line. The geometrical properties 
of the line as a function of length-scale can be conveniently described using 
the instantaneous quadratic width: 



w'it) = jj:[uix,t)-uit)]\ (5) 

-^ x=0 



or the averaged structure factor 

(6) 



L-l 2 






5 



where q = 2Tm/L, with n = 1, . . . , L — 1, and u{t) = L~^ Ylx=o ^(^' ^) i^ the 
center of mass position. When the steady state regime is reached, these quan- 
tities become time independent. In particular, for the three reference states 
(equilibrium, depinning and fast flow), the last two quantities display self-affine 
behavior: w'^ ~ L^'' and Sq ~ l/g"'"'"^^. For different values of the external force, 
the above observables display crossovers between the corresponding reference 
states. 



4 Numerical methods and main results 



At zero force the steady state is the equilibrium state which corresponds to 
the ground state at T = 0. There exist exact and efficient numerical meth- 
ods for targeting the equilibrium both at zero and finite temperature. These 



3. For developing some algorithms and specially for the statics, it might be also 
convenient to discretize the displacement field u{x,t). 



methods rely on transfer matrix techniques and apply to short range elas- 
ticity [37]. Upon increasing the driving force from zero, we enter first in the 
creep regime which is defined at low, but non zero temperature. In this regime 
the Langevin dynamics can be implemented at finite temperature [38,39], but 
becomes highly inefficient in the vanishing temperature limit (T — )> 0^). For- 
tunately, in this limit an exact algorithm allows to characterize the ultra slow 
motion of the line [31,32]. For a larger value of the external force we reach the 
depinning transition (/ = /c). At T = direct Langevin dynamics simulations 
converge very slowly to the correct steady-state. Different methods have been 
thus developed: (i) an efficient algorithm allows an exact calculation of the 
critical force and critical configuration for each disorder realization [40]. (ii) 
Molecular dynamics simulations of the non-steady relaxation allow to get the 
dynamical critical exponents and target the thermodynamic critical force [41]. 
At T > and / = /c, Langevin dynamics allows to describe the critical 
thermal rounding of the depinning transition [42,43]. Finally simple Langevin 
dynamics in the large force limit (/ ^ /c) captures the features of the fast 
flow regime at all temperatures. 



4-1 Equilibrium State 



At zero force, the geometrical properties of the elastic string can be targeted 
using a transfer matrix method for the discrete directed polymer model. The 
line is described by the discrete variable u{x), that gives the displacement of 
the string on the slice x of a square lattice, and a disorder potential U, which 
is drawn on each site (x, u). The energy of a given configuration is given by the 
sum of the site energies along the path u{x): E = J2x U{x, u{x)). Furthermore, 
the so-called solid-on-solid (SOS) restriction \u{x + 1) — u{x)\ = 1 has to be 
implemented. When / = this model belongs to the same universality class 
as Eq. (1). 

Given a disorder realization, the ground state configuration of a polymer start- 
ing in (0, 0) and ending in (x, u) can be found using the following recursive 
relation: 

^min(a;, u) = U{x, u) + min[E^i„(x - 1, m - 1/2), E^^^ix -l,u + 1/2)]. (7) 

At finite temperature T it is possible to compute the weight Z^^u of all polymers 
starting in (0,0) and ending in {x,u) using the following recursion: 



Zx,u — e ^'"^ \Zx-l^u-l/2 + Zx-l^u+l/2J 



with the initial condition Zo,„ = (5o,u. Therefore, the probability to observe 
a polymer ending in (x,m) is ^x,u/Z]«' ^x,u'- Due to the recursion relation. 
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Figure 2. Low temperature dynamics of the driven elastic string below the depin- 
ning threshold. The optimal path to escape from a given metastable configuration 
a pass through a saddle configuration /3, from which the system is able to relax 
deterministically to the next metastable configuration 7 with E^/ < E^. 

Zx^u grows exponentially with x. To avoid numerical instability, all weights 
Z^^u at fixed x have to be divided by the largest one, which does not change 
the polymer ending probability [44]. With these procedures it is possible to 
characterize both the geometrical properties of the polymer (as the rough- 
ness exponent Ceq = 2/3 in Id) and the free energy or ground state energy 
fluctuations (as the characteristic exponent 9) [45,46]. 



4-2 Creep regime 



For non-zero drive but still below the depinning threshold (0 < / < /c), a 
finite temperature must be imposed in order to forget the initial condition 
and reach a unique stationary state. In references [31,32] the regime relevant 
for creep motion was studied; this is, the steady state regime of an interface 
in the limit of vanishing temperature. A direct numerical study of the motion 
of such regime is very difficult. The motion takes place by thermal activation 
over barriers leading to extremely long activation times making numerical 
techniques such as the Langevin dynamics inefficient. Fortunately, different 
numerical method [31,32] allows to follow the motion of an interface at finite 
temperature, avoiding the above-mentioned difficulty. This method directly 
implements the polymer dynamics as a unique forward-moving sequence of 
metastable states of decreasing energy. It can be proved that such a sequence 
corresponds to the exact dynamics for a finite interface in a given disorder 
realization in the T — )■ limit. 



The elementary step of such a dynamics is sketched in Fig. 2. It consists in 
finding the optimal path from one metastable configuration a with energy 
Ea to the next metastable configuration 7, with a lower energy. This path 
pass through a saddle configuration [3 with an energy Ep > E^- Once in /3 
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Figure 3. (a) Steady state structure factor S{q) of the line in the T — t- hmit for 
different forces (curves are shifted for clarity), (b) The steady state properties of 
the elastic string at T ^- are determined by £opt (filled symbols) and by £av These 
lengths separate regions characterized by the equilibrium exponent, the depinning 
exponent (gray region), and the fast flow one. The divergent length ^reiax (open 
symbols) is associated only with transient dynamics. Lines are guides to the eye. 
Figures are adapted with permission from [31]. 

the polymer relaxes deterministically to 7. The configurations a and 7 differ 
from each other by a portion of size ^reiax, while a and P differ by £opt- The 
maximun energy E^^^ reached by the polymer through the path defines a drive 
dependent barrier Eb{f) = -Emax — E^, therefore the time associated to this 
elementary step is given by ~ e^''^^^^'^ . 

Numerical results are summarized in Fig. 3: (i) The structure factor behavior 
below /c supports the picture proposed in Fig. 1. In particular this behavior 
is consistent with FRG calculations [16] that predict in the creep regime the 
existence of a characteristic scale (namely £opt) below which the system is 
at equilibrium and above which the system is characterized by deterministic 
forward motion. Moreover, in qualitative agreement with these predictions it is 
observed that both i^pt and Eh{f) increase as / decreases. Unfortunately this 
algorithm has not allowed to explore the region of vanishing / where iopt is 
expected to diverge as ~ y-i/(2-Ceq) g^nd the "optimal" barrier as Eh(f) ~ f~^. 
In turn, it can be seen that the size irei diverges as ~ (/c — /)~'^'*''p and this 
is interpreted as the linear size of the deterministic depinning-like avalanche 
triggered by thermal nucleation. It is worth remarking that although freiax 
diverges as / — )■ /c from below, it does not affect the steady-state spatial 
correlations, unlike £av and iopt- 

In order to study the effects of larger temperatures one can use Langevin 
dynamics [39]. These simulations show that £av has a finite value for T > 
and < / < /c and it actually diverges as / — )■ (as sketched in Fig. 1), but 
the lack of precison does not allow to test the universal behavior £av ~ j"-" j-^^ 
predicted by FRG calculations [16]. Note that £av diverges as the velocity goes 
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Figure 4. Results for the critical state in finite samples obtained with the exact 
algorithm of Sec. 4.3.1. (a) Critical configuration structure factor (see definition 
in Eq. (6)). The fitted roughness exponent yields Cdep = 1.250 it 0.005. (b) The 
finite-size critical force /c as a function of the longitudinal size L for periodic samples 
of size L X M with M = kL^'^'^f . As L grows the asymptotic value of the critical 
threshold becomes independent on the aspect ratio k. The dashed line corresponds to 
the thermodynamic limit fc = 1.916 it 0.001. Figures are reproduced with permission 
from [41]. 

to zero, while freiax remains finite in the same hmit. In other words, while in 
the T — )■ limit the string is blocked in the first metastable state with lower 
energy producing an avalanche of typical size ^reiax, at finite temperature a 
larger avalanche of size £av takes place. In the creep regime these depinning- 
like processes are fast in comparison with the waiting times for the thermally 
activated jumps (^opt)- 

Finally, the finite temperature long time relaxation of a flat line in absence 
of drive can be described as a creep process over barriers that increase in 
time, such that the growing correlation length (separating the equilibrated 
length-scales with the ones retaining memory of the initial condition) is [47] 
£(t)~[Tlog(t)]V^ 



4-3 Depinning regime at T = 



4-3.1 Critical force and critical configuration 

At T = 0, the transition between a pinned phase and a moving phase is 
displayed in any finite sample with a given disorder realization. This is assured 
by a set of properties pointed out by Middleton [48] . The first one is the so- 
called "no-passing rule": if two strings u{x,t) and u{x,t) do not cross at a 
given time, they will not cross at any later times. Another important property 
of Middleton's theorem states that if, at an initial time, the velocities are 
non-negative for all points x, they will remain non-negative for all later times. 
It follows from these properties that, once we have found a forward moving 
string, we can be sure that snapshots of the string at later times will be ahead 
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(in the u direction). As a consequence of these two properties, we can define 
(for a given finite sample) the critical force f^^P^^ as the maximal force for 
which a metastable configuration still exists; i.e., a configuration for which al 
points X verify dtu{x, t) = 0. Analogously to the ground state polymer for the 
statics (Sec. 4.1), this last unique metastable configuration displays self-affine 
properties with a well defined roughness exponent Qep (see Fig. 4a). For a 
finite sample, f^^P^^ is a stochastic variable which depends on the disorder 
realization. As the size L of the sample grows, the asymptotic critical threshold 
fc is approached (see Fig. 4b), and the fiuctuations of J^^p^^ are suppressed 
as 



/■sample _ /.sample ^ j^-2/uiep. (Q\ 



where z/^ep is the characteristic correlation length exponent. Both the expo- 
nents (dep and Udep can be carefully determined using numerical algorithms 
that avoid the direct numerical integration of Eq. (1). The key idea behind 
these methods is to check for the existence of metastable states for the string 
in the particular sampled disorder for a given value of the external force /. 
This check can be performed in a computing time that grows only linearly 
with the system size using properties discussed by Middleton [40,50]. 



4-3.2 Critical and dynamical exponents: the non-steady relaxation 

The accurate determination of critical exponents is a difficult task, specially for 
critical phenomena in disordered systems. The so-called Short-Time Dynamics 
method (STD) allows to obtain critical exponents values without the need of 
system equilibration. It relies on the validity of an homogeneous relation for 
the order parameter when the system performs a non-steady relaxation at the 
critical point. This scaling relation, in contrast with traditional equilibrium 
finite-size scaling methods, can be used to avoid finite-size effects but includes 
dependences with time and initial conditions [51,52,53]. 

To apply the STD to the depinning transition we rely on the analogy with 
standard phase transitions [36]. By considering the velocity v as an order pa- 
rameter, and the dimensionless force {f — fc)/fc as the reduced driving field, we 
expect the following homogeneity relation to be valid for the long-time relax- 
ation of an initially fiat string (infinite velocity initial conditional) [54,55,41]: 

v{h, L, t) = t-^/'''d±{t^/"'h, t-^/'L); (10) 

where h = \f — fd/fc and the function v± has two branches depending on 



4. This is equivalent to the case of an "ordered" initial condition for the order 
parameter. 
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7 1.4 




Figure 5. Evolution of 7(t) (a) and v{t) (b) for / « fc, with fc the thermody- 
namic critical force. Two other forces, just above (f^y) and just below (/~) fc are 
also shown. Continuous lines are fits for 7(t) and v{t) including power-law correc- 
tions to the asymptotic scaling yielding the true critical exponents. In (a) vertical 
dashed lines separate qualitatively different time regimes. The inset in (b) shows 
that corrections to scaling are power-law like. Figures are adapted with permission 
from [41]. 

the sign of / — fc- Close to depinning, the relaxational dynamics described 
by Eq. (10) is valid in a short-time regime, after which the velocity reaches 
a steady-state value if / > fc, while for f < fc the string is blocked in a 
metastable state (with memory of the initial condition) and v ^ 0. Exactly 
at the critical point {h = 0) and in the limit L — )■ oo we expect a power- 
law behavior for the velocity, v ~ t~l^/'^^. Alternatively, we can write v ~ 
(,{t)~^/'^ ., where i,{t) ~ t^^^ is the dynamical length, describing the growth of 
roughness as w{t) ~ (i{t)^ ~ t^^^ . By fitting v{t) and w(t) we thus in principle 
get P/vz and Q/ z independently. Since C (and therefore v) can be obtained 
independently using the method described in Sec. 4.3.1, with the STD method 
it is possible to extract in addition the critical exponent /3 and the dynamical 
exponent z. 

In Ref. [41] large scale Langevin dynamics simulations have been implemented 
to obtain v{t) and w{t) over a large period of time for large system sizes, 
such that L ^ i{t). This has allowed to detect robust (under different shape 
and nature of the disorder correlator) scaling corrections to the asymptotic 
power-law forms, and thus to avoid fitting effective power-law decays, yield- 
ing significantly biased incorrect exponents. A practical criterion to estimate 
the crossover between the "mesoscopic" time regime with corrections and the 
truly universal "macroscopic" time regime was found by observing that the 
hyperscaling relation /3 = i>{z — Q is violated in the former regime. This can 
be easily seen by analyzing the quantity 7(t) = ^|4: if the measured 7(t) 
depends on time, it means that the system has not reached the critical scaling 
regime yet. 

In Fig. 5(a) we show the behavior of 7(t). After a microscopic regime, 7(t) 
slowly decreases implying an "effective" unbalanced relation [u^z — C) ~ /^]es ^ 
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Figure 6. (a) Structure factor Sg at a finite temperature T = 0.02 showing two 
regimes with different roughness exponents, Cff and Cdep, above and below the 
crossover length scale ^av, respectively, (b) Velocity- force scaling function in the 
thermal rounding region. The disorder intensity is kept fixed and the velocity for 
different temperatures collapses in a single function. Figures are adapted with per- 
mission from [43] . 

0. In this mesoscopic regime the critical dynamics scaling is not valid. After a 
surprisingly long crossover however, 7(t) develops a plateau. This information 
allows to determine the appropriate time regime of v{t) and w(t) where it is 
possible to fit (3/h'z and (/z, respectively. In Fig. 5(b) the behavior of v{t) is 



1 + [t/t,] 



(with 



shown. Power- law corrections of the form v{t) = Vq t '^'^^ 
Vo, ty, and a^ fitting parameters) are find. At large times, when 7(t) develops a 
plateau, an accurate power-law fit can be performed for v{t) at the thermody- 
namic depinning force fc = 1.5652 ± 0.0003, yielding 13 /{vz) = 0.128 ± 0.003. 
The observation of scaling corrections explain the systematically larger effec- 
tive values for (P/uz) reported before for smaller systems [41]. An analogous 
analysis for w{t) shows similar robust power- law scaling corrections, and gives 
a consistent value for the exponent, 1 — (/z = 0.128 ± 0.003. These results, 
combined with an independently determined value oi ( = 1.250 ± 0.005 (see 
Fig. 4a) and the relation u = 1/(2 - (), allows to get /3 = 0.245 ± 0.006, 
z = 1.433 ±0.007, C = 1.250 ±0.005 and z/ = 1.333 ±0.007 [41]. Scaling forms 
for the joint force-time dependence of the velocity extracted from Eq. (10) are 
in excellent agreement with these exponents [41]. 



4-4 Depinning regime atT>0: the critical thermal rounding 



When the temperature is finite, there is no sharp transition between zero and 
finite velocity regimes. At forces around the critical value, / ~ /c, a finite 
temperature value smears out the transition, which is no longer abrupt. This 
thermal rounding of the depinning transition can be characterized, exactly at 
the critical force / = fc, by a power-law vanishing of the velocity with the 
temperature v ~ T'^, where ip is the so-called the thermal rounding expo- 
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nent [38,42,43,56,57]. 

At small length scales, q ^ 1/^av, the structure factor shows the typical rough- 
ness regime associated to depinning, i.e. Sq ~ g-(i+2Cdep)^ while at large length 
scales, q ^ l/^av, geometrical properties are dictated by effective thermal fluc- 
tuations induced by the disorder, i.e. Sq ~ g-{i+2Cff)^ ^s shown in Fig. 6(a). 
In the critical region the depinning correlation length is given by the velocity 
as £av ~ v~^/^ . Thus, the depinning correlation length depends on the tem- 
perature only through the velocity and in the thermal rounding regime [42] 
£av ~ T"^*^/^. With this information one can write for the structure factor 



where the scaling function s{y) ~ y-{i+2Cff) for ?/ ^ 1 and s{y) ~ y-{i+2Cdep) 
for y ^ 1. In Ref. [42] it was shown that the structure factor scales with the 
previous form. 

All the information about the thermal rounding regime can be gathered in the 
expected universal behavior oiv{h, T). It is possible to show that assuming the 
homogeneous relation between the velocity and both temperature and force, 
as it is usual for phase transitions, a universal function follows. If there were 
not strong finite size effects, in the vicinity of the critical region the velocity 
should scale as 

vT-^ ~ h± (/T-^/'^) , (12) 



with h±{x) a universal function which can in principle be different above (/i_|_) 
and below (h-) threshold. One expects that h+ ~ y^ for ?/ ^ 1. Figure 6(b) 
shows this scaling relation for fixed disorder intensity and different temper- 
atures as indicated. This universal behavior has been numerically confirmed 
using Langevin dynamics simulations of the QEW equation [43] and indirectly 
tested in experiments of domain wall motion in 2D ferromagnets [58] . 



5 Conclusions and Perspectives 



The driven QEW model represents a minimal paradigmatic model for studying 
the dynamics of driven elastic systems in random media. In this paper we 
have reviewed a series of numerical methods specially developed for studying 
different and key properties of this system. The results found in the literature, 
lead us to the physical picture discussed in Sec. 2 and ilustrated in Fig. 1. 

The most challenging aspect of the problem remains the quantitative compar- 
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ison with experiments. On one hand, it would be important to make quanti- 
tative predictions for different elastic universality classes. Although we expect 
the same basic physics described here to emerge at large enough scales, it 
is important to accurately obtain the corresponding critical exponents and to 
understand the additional dynamical crossovers that may arise at intermediate 
length-scales. 

At last, probably the least understood but still experimentally relevant phe- 
nomena for the dynamics of interfaces in random media are: the occurrence 
of "plasticity" (displayed by overhangs and bubbles [59]), the effect of inter- 
nal degrees of freedom (such as the "spin phase" coupled to the position of 
magnetic domain walls [5,60]), and the effect of "structural relaxation" [15] 
in host materials. This may encourage the development of new models and 
novel efficient numerical approaches in the field. 
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